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************************* Least Square Lattice ********«•»»»«•••«*********** 
************************ Noise Cancelling ****************************/ 

/* Example for constant saturation approach to noise cancelling */ 
^define LAMBDA 0.95 

void OxiLSL_NC( int reset, 

int passes, 

int sat_f actor, 

int *signal_l, 

int *6ignal_2, 

int *target_l, 

int *target_2) { 

int i, ii, k, m, n, contraction; 
static int *s_a, *s_b, *out_a, *out_b; 
static float Delta_sqr, scale, noise_ref; 

if ( reset *~ TRUE) { 
s_a *= signal_l; 
s_b « signal_2; 
out_a « target_l; 
IM= out_b = target_2; 
m scale « 1.0 /4160.0; t - 

y /* noise canceller initialization at time t=0 */ 

+= 

s nc[0].berr = 0.0; 

Q nc[0] .Gamma «= 1.0; 

~J for(m*0; m<NC_CELLS; m++) { 

H nc[m].err_a = 0.0; 

N nc[mj.err~b « 0.0; 

y3 nc[m].Roh_a = 0.0; 

gi nc[m].Roh_b = 0.0; 

~ nc(m). Delta = 0.0; 

nc[m] .Fswsqr = 0. 00001 ; 

nc(m].Bsvsqr = 0.00001; 



************************************** 



} 



END INITIALIZATION ■■' ■ «—«««*/ 



for(k«=0; k<passes; k++) { 



contraction - FALSE; 

for(m»0; m< NC_CELLS; m++) { /* Update delay elements 

nc[m].berrl « nc[m).berr; 
nc(m] .Bswsqrl - nc [m] . Bswsqr ; 

} 

noise ref - sat_f actor * log (1.0 - (*s_a) * scale) 

- log (1.0 - (*s_b) * scale) ; 
nc[0],err a « log(1.0 - (*s_a) * scale); 
nc[0].err~b « log(1.0 - (*s_b) * scale); 

++s a; 



♦ 




a ; 



i s 



++s_b; 

nc[0] . ferr 
nc[0] .berr 
nc[0] .Fsvsqr 
nc[0] .Bswsqr 



noise_ref ; 
noise_ref ; 

LAMBDA * nc[0]. Fswsqr + noise_ref * noise_ref; 
nc[0] .Fswsqr; 



/* Order Update */ 

for(n«=l;( n < NC_CELLS) ££ (contraction 

/* Adaptive Lattice Section */ 

m - n-1; 
ii- n-l; 



FALSE) ; n++) { 



nc[m]. Delta *■ 
nc[mj .Delta 
Delta_sqr 

nc[n].fref 
ncfnj.bref 



LAMBDA; 

nc[m).berrl * nc[m].ferr / nc[m] .Gamma ; 

nc[m). Delta * ncfm]. Delta; 

-nc[m). Delta / nc[m] .Bswsqrl ; 

-nc(mj. Delta / nc[m] .Fswsqr; 



nc[n).ferr = nc[m).ferr + nc[n).fref * nc[m].berrl; 
nc(nj.berr = nc[m].berrl + ncfnj.bref * nc[m].ferr; 

nc[n). Fswsqr = nc[m]. Fswsqr - Delta_sqr / nc(m) .Bswsqr 1; 
nc[n). Bswsqr = nc [m] .Bswsqrl - Delta_sqr / nc[m] .Fswsqr; 

if( (nc[n] .Fswsqr + nc[h] .Bswsqr) > 0.00001 | | (n < 5) ) { 

nc(n]. Gamma «= nc{m]. Gamma - nc[m).berrl * nc[m).berrl / nc[m] .Bswsqrl 
if (nc[n] .Gamma < 0.05) nc[nj. Gamma «= 0.05; 
if (nc[nj .Gamma > 1.00) nc[n]. Gamma « 1.00; 

/* Joint Process Estimation Section */ 

nc(m].Roh_a LAMBDA; 

nc[m].Roh_a += nc{m].berr * nc[m].err_a / nc[ra]. Gamma ; 
nc(m].k_a = nc[m].Roh_a / nc[m) .Bswsqr ; 
nc(n).err_a = nc[m].err_a - nc[m].k_a * nc(m].berr; 

ncfm).Roh_b *= LAMBDA; 

nc[m].Roh_b += nc[m).berr * nc(m].err_b / nc[m]. Gamma ; 

nc[m).k_b = nc[m).Roh_b / nc[m] .Bswsqr ; 

nc[n].err_b = nc[m].err_b - nc[m].k_b * nc[m].berr; 

} 

else { 

contraction «= TRUE; 
for(i«=n; i<NC_CELLS; i++) { 



nc[i] . err_a 
nc[i] .Roh_a 
ncfi] .err_b 
nc[i] .Roh_b 
nc[i] .Delta 
nc( i ) . Fswsqr 
nc[i] .Bswsqr 



0.0; 
0.0; 
0.0; 
0.0; 
0.0; 

0. 00001; 
0.00001; 



nc(i] .Bswsqrl = 0.00001; 



) 

) 

t 

) 

*out_a++ = (int)( (-exp(nc[ii).err_a) +1.0) / scale) ; 
*out_c++ » (int) ( (-exp(nc[ii].err_b) +1.0) / scale) ; 

} 

) 

/ft********************* Least Square Lattice 



